Spatial structures and dynamics of kinetically constrained models of glasses 
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Kob and Andersen's simple lattice models for the dynamics of structural glasses are analyzed. 
Although the particles have only hard core interactions, the imposed constraint that they cannot 
move if surrounded by too many others causes slow dynamics. On Bethe lattices a dynamical 
transition to a partially frozen phase occurs. In finite dimensions there exist rare mobile elements 
that destroy the transition. At low vacancy density, v, the spacing, 3, between mobile elements 
diverges exponentially or faster in 1/v. Within the mobile elements, the dynamics is intrinsically 
cooperative and the characteristic time scale diverges faster than any power of 1/v (although slower 
than 3). The tagged-particle diffusion coefficient vanishes roughly as H _d . 



In many liquids dramatic slowing down occurs on cool- 
ing, equilibrium cannot be achieved, and structural relax- 
ation becomes complicated and spatially heterogeneous 
Despite a great deal of effort, these remarkable 
phenomena associated with the "glass transition" are still 
far from understood. Indeed, the most basic issues are 
unresolved: Is the rapid slowing down due to proximity 
to an equilibrium phase transition, albeit in a restricted 
part of phase space? Or is the underlying cause entirely 
dynamical? In cither case, are there characteristic length 
scales that grow substantially near the glass transition? 
Theoretical developments have been hampered by a 
shortage of models that capture essential features yet are 
simple enough to be analyzed. Exceptions are kinetically 
constrained models 0, based on the Ansatz that the 
glass transition is caused by geometrical constraints on 
dynamical rearrangements with static correlations play- 
ing no role (see e.g. 6]). These are also models for gran- 
ular media (see e.g. 7]) for which slow dynamics occurs 
already at densities well below close-packing 
We focus on one of the simplest kinetically constrained 
models, the Kob- Andersen model (KA) 9] . KA is a single 
component lattice gas with no static interactions other 
than hard core exclusion and dynamics given by a contin- 
uous time stochastic process. Each particle attempts, at 
a fixed rate, to move to a randomly chosen empty neigh- 
boring site; but the jump is allowed only if both before 
and after the move the particle has no more than some 
number, to, neighboring particles. This corresponds to 
vacancies moving only if the initial and final sites have 
at least s = z — to — 1 neighboring vacancies, with z the 
coordination number of the lattice. Since this dynamics 
satisfies detailed balance, the trivial distribution that is 
uniform over all configurations with a fixed number of 
particles is stationary: there can thus be no equilibrium 
transition. Nevertheless, the dynamics is sluggish and 
heterogeneous at high density |9j, [l0| . For a three dimen- 



sional cubic-lattice with s — 2, fits of the self diffusion 
coefficient of a tagged-particle, Ds, strongly suggest a 
dynamical glass transition at p c ~ 0.881, above which 
Ds appears to vanish and the structural relaxation time 
to diverge 

Our analysis shows that there are four classes of behav- 
ior for general KA models. These can be characterized 
by the dependence of Ds on the vacancy density, v: (N) 
normal with Ds ~ v q as v — ► 0; (C) collective freezing 
with Ds — > faster than any power of v as v — > ; 
(T) dynamical transition with Ds — > as v \ V c with a 
non-zero critical v c — 1 — p c ; and (F) always frozen with 
a finite fraction of the particles never moving for any v. 
Normal behavior occurs if a finite cluster of q vacancies 
can move through an otherwise-totally filled lattice, e.g. 
a triangular lattice with s = 1 in which vacancy pairs 
(q = 2) are mobile. This yields Ds proportional to the 
cluster concentration, v q . At the opposite extreme, a 
square lattice with s = 2 is always frozen as a four par- 
ticle square can never move. We focus on the interesting 
intermediate cases, C and T. We show that a dynamical 
transition takes place on (tree-like) Bethe lattices — even 
with finite size loops. But such a transition cannot oc- 
cur on finite dimensional lattices: if neither normal nor 
always frozen, these exhibit collective freezing with the 
dominant dynamical processes involving a number of va- 
cancies that diverges as v — ► 0. This class (C) includes 
the original s = 2 cubic lattice. 

We first analyze Bethe lattices: infinite tree-like graphs 
with fixed connectivity z which crudely approximate 
high-dimensional or high-coordination-number lattices. 
Monte-Carlo simulations of 10 sites with z = 4 and s = 1 
(crudely mimicking the s = 1 square lattice) suggest a 
dynamical correlation length (as defined i n 1121 1 that di- 
verges at a critical density p c ~ 0.89. 11] Approach- 
ing this apparent dynamical transition the local-density 
correlation function displays two-step relaxation, as for 



supercooled liquids. 

The tree-like structure of Bethe lattices enables ana- 
lytic study by iteration. Arranging the tree with k = z— 1 
branches going up from each node and one going down, 
we focus on the following events for a chosen node and 
rearrangements restricted to the subtree above it : (A) 
node is occupied by a particle which cannot move up 
even if the site below it is empty; (B) node is empty but 
a particle can never move onto it from below; (C) node 
is occupied by a particle which can only move up if the 
node below it is empty; (T>) none of the above. Self- 
consistent equations for the probabilities of these events 
can be found by combining k branches together with the 
site below them to which they are all linked For 
,s = 1 with z = k + 1 = 4, this yields a critical den- 
sity, p c = 0.888825 in agreement with simulations. Be- 
low p c , Prob[_4] = Prob[S] = 0; particles are mobile 
and diffusive at long times. At p c , Prob[.4] = A c ^ 
and the fraction of immobile particles jumps discontin- 
uously; so does the extensive configurational entropy, 
Sc = ln[#different ergodic components]. 0] Above 
Pc Prob[^4] — A c r~j y ' p — p c , the singularity indicative 
of diverging length and time scales. The KA dynam- 
ical transition on this Bethe lattice (like 1/r 2 percola- 
tion in one-dimension [l^l thus has both critical and 
"first order" characteristics. Interestingly, infinite-range 
quenched random p-spin models, conjectured to be re- 
lated to the glass transition, also exhibit a dynamical 
transition with similar features 

The basic features found above hold for generic KA 
models on Bethe lattices with s ^ 0, k, including those 
with finite size loops [ll|> i- e - f° r the so called pure 
Husimi trees and for Bethe lattices where single sites 
are replaced by L x L squares (these interpolate among 
the Bethe lattice, L = 1, and the square lattice, L = oo). 
The only exceptions are cases in which particles with two 
occupied neighbors cannot move (to =1): the transition 
is then continuous and p c = l/k, the critical density for 
conventional site percolation. This equivalence is due 
to the fact that the dynamical transition is associated 
with the emergence of an infinite cluster of permanently 
blocked particles. In general a frozen phase will only exist 
if an infinite cluster remains after all the particles mobile 
under KA rules have been iteratively removed. This is a 
type of bootstrap vercolation, already introduced in this 
[tj and broader llalrRll8| contexts. For Bethe lattices, 
an exact solution |ll| fsee |li| for a simpler example) 
yields the same critical density found above. 

We now turn to finite dimensions and show that the dy- 
namical transition is destroyed by exponentially rare pro- 
cesses which will only occur in sufficiently large systems. 
Nevertheless, at high densities the dynamics is intrinsi- 
cally collective and extremely slow. Furthermore, the 
ghost of the Bethe lattice transition may cause a sharp 
crossover at an apparent critical density. 

We focus initially on the simple s = 1 square lattice 



case. First, let us define a configuration as framed if all 
its boundary sites are empty (see for example the square 
on the right of fig. 1); and as frameable if, by an al- 
lowed sequence of moves, a framed configuration can be 
reached (see for example the square on the left of fig. 1). 
A key observation is that the lines of vacancies on the 
boundaries of a framed square can be shifted (moving 
vacancies starting from the corner) so that any nearest 
neighbor pair of particles inside the square can be sand- 
wiched between two lines of vacancies. The position of 
this pair of particles can then be exchanged and the va- 
cancy lines afterwards shifted back to their original posi- 
tions. The net result is a pair exchange of particles with 
all other particles returned to their original positions. By 
combining such pair exchanges all frameable configura- 
tions with the same density can be connected. Thus, 
these form an irreducible component Z p of the configu- 
ration space. We must now show that almost all random 
configurations with a given p < 1 belong to Z p in the 
thermodynamic limit. This can be done iteratively by 
observing that an I by I frameable configuration that has 
at least two vacancies externally adjacent to each of its 
sides is also an I + 2 by I + 2 frameable configuration, 
see fig. 1. Starting from a two by two frameable "nu- 
cleus" of vacancies, one can thus grow an L by L frame- 
able configuration if the requisite vacancies are present 
in each concentric shell. The probability that this occurs 

is p L (p) - (i - P f nf=4 a [i - p 21 - 2/p 2/ - i (i - p)] 4 which 

converges to a non-zero probability, Poo (p). Thus for all p 
infinite frameable squares exist with all frameable config- 
urations within them being reachable. For small vacancy 
densities, Poo(p) ^ e~ 2KF/v with K F ~ 4.48 . Moreover, 
Pl(p) depends weakly on L for L > £(p) ~ ln(l/v)/v; £ 
is thus the core size of frameable regions. Note that Pl 
is the probability that a frameable square can be con- 
structed around a nucleus at a fixed position. On the 
other hand the probability, Fl(p), that a large square is 
frameable is roughly given by considering the <3([L/£] 2 ) 
possible positions of a nucleus within the square. Al- 
though the probabilities of a configuration being frame- 
able around each of these possible nuclei are not inde- 
pendent, if the nuclei are separated by more than a core 
diameter they are roughly independent. Indeed, it can 
be proven that Fl asymptotes to one for L ^> 5 with the 
long crossover length scale S(/?) given, as expected from 
the above argument, by E 2 P^(p) ~ E 2 P OQ (p) ~ £ 2 . 

We have shown that at any density almost all con- 
figurations of sufficiently large systems can be reached 
from one another. This excludes the possibility of a dy- 
namical transition [2(J • The phase space of systems with 
linear size, L, larger than the crossover length 3 is cov- 
ered almost entirely by a single ergodic component, while 
the phase space for L < S is typically decomposed into 
many disjoint parts. The framing analysis yields an upper 
bound for S. 
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FIG. 1: On the left there is a 4 by 4 framed configuration 
with two vacancies (white circles) adjacent to each side. Af- 
ter moving the external vacancies at the corners the framed 
square can be expanded to the 6 by 6 framed configuration 
at right. 



A lower bound for S can be obtained from bootstrap 
percolation arguments. ^From a random configuration 
with density p iteratively remove all particles with less 
than three occupied neighbors. If any particles remain 
they must form a system-spanning cluster of frozen par- 
ticles and the configuration space will thus be broken up 
into many disjoint pieces. This happens with high proba- 
bility for InL < Kp/v with K F = tt 2 /18 003, yielding 
a lower bound for S with the same density dependence as 
the above upper bound, but a different constant Kf less 
than Kf of the upper bound. However, a better defini- 
tion of framed configurations yields a Kf which is identi- 
cal to Kp jll) . The corresponding frame is minimal and 
the asymptotic form of the crossover length S ~ e n / 18u 
exact (up to subdominant factors). 

All other finite-dimensional KA models can be ana- 
lyzed by generalizing the notion of frameable configura- 
tions, we thereby rule out any dynamical transitions. For 
rf-dimensional hyper-cubic lattices s = and s > d — 1 
correspond respectively to N and F while all the non- 
trivial cases, 1 < s < d—1, are class C, collective freezing. 
Analogous expansion arguments yield upper bounds 

on their crossover lengths: H < exp os (Ku(s, d)/v' J ^^j 
with exp os the exponential function iterated s times and 
Ku(s, d) a positive constant given in 1 1 11] - Bootstrap per- 
colation yields a lower bound of the same form but with a 
different constant Kl(s 1 d) [2l| The scale, £, above which 
the probability that a nucleus is expandable saturates, is 
generally of order In S, up to log-log factors. 

We now turn to the dynamics, focusing again on the 
s = 1 square lattice model at high density. Vacan- 
cies will typically be far apart or in small clusters that 
cannot move. However vacancies can move within the 
cores of large frameable regions. And these cores of size 
£ ~ ln(l/v)/v are themselves mobile because they are 
likely to find vacancies on all the £ successive line seg- 
ments needed for them to move in a given direction. A 
generic particle cannot move substantial distances ex- 
cept when a mobile core passes by; assuming indepen- 
dent core motion a tagged particle will thus diffuse with 



Dg ~ Dj^um, where nu ~ 1/S 2 is the density of the 
mobile cores which diffuse with Dm — £ 2 / t ?j V T £ being 
the typical relaxation rate of a core. Using our framing 
analysis and a generalization of the technique of , we 
have proved that, indeed, D$ > for all p. 

The relaxation time, tj, of a core is dominated by 
the time to reach the most severe bottleneck in config- 
uration space. Since all configurations are equiprobable 
this is proportional to exp(AS'), the ratio of the num- 
ber of accessible (frameable) configurations of the core 
to the number in the bottleneck [llj. The worst case 
scenario in which there is a single bottleneck configu- 
ration corresponds to AS = Stotai — ln£!, yielding an 
upper bound T£ < £!. A lower bound is obtained by 
noting [n] that to equilibrate a frameable square of size 
£ one has to pass through configurations with the nu- 
cleus in a corner of the square. Using a transfer ma- 
trix technique we obtain for £ x £ minimally frameable 
squares an entropy difference AS ~ T\/I + a\n£ + C 
with T = 2v / 6~+V 7 22 - 2^3 = 3.075 and a computable 
in principle. Simulations of frameable squares with 
I — 4 to £ — 16 yield equilibration times with log re in- 
creasing slower than I and in good agreement with the 
square root law with the predicted coefficient T plus log- 
arithmic corrections. We thus conjecture that the en- 
tropic bottlenecks for moving the core of a mobile region 
are equivalent, up to subdominant factors, to the entropy 
loss associated with forcing the nucleus of of such core to 
a corner. 

As long as bottlenecks for motion of the mobile cores 
are not too much larger than the above estimates, T£ <C 
S 2 and the dominant contribution to D$ ~ DmTim will 
arise from the low density of mobile regions rather than 
the long time needed for motion within them. Thus, 
to leading order in In Ds ~ 2Kf/v for the square 
lattice model with s = 1. Numerical simulations at 
high densities fit this form well, see fig. 2: in particu- 
lar lim^o^hi-D = 2Kp — 0.9 — 0.95 cf the predicted 
2Kp = 1.1. For the other KA models, we obtain similar 
relationship between Dg and Hm- For the s = 2 cubic 
case theis yelds lnlnl?^ 1 oc X/v , a result consistent with 
a recent finite size scaling analysis [23). 

Thus far we have focused on the high density limit. 
Yet simulations on KA and other kinetically constrained 
models suggest a transition at a non-trivial critical den- 
sity 0- This is a natural consequence of the existence 
of a dynamical transition on Bethe lattices: in finite di- 
mensions, the transition will be replaced by a crossover 
which, in some cases, could be quite sharp. In the s = 1 
square lattice case a three-vacancy element can move 
along a network of other vacancies provided these are 
linked no more weakly than via second neighbors (along 
axes or diagonals). At low density, there will always be 
a percolating vacancy network of such type. As the den- 
sity increases, diffusion slows down as the vacancy per- 
colation cluster, on which mobile elements rely, sparsens. 
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FIG. 2: (InD)" 1 as a function of p for 400 by 400 square 
lattice with s = 1. Straight line shows In Ds oc 1/v behavior. 

Beyond the second-neighbor vacancy percolation transi- 
tion at P2P, the simple mobile elements can no longer 
move. But by growing, mobile elements can still form: 
they have to become large enough to find a percolating 
cluster on which objects of their size can move. At high 
densities, such mobile elements will be the cores of (min- 
imal) frameable regions that we have discussed above. 
If the crossover from small mobile objects and behavior 
similar to that on Bethe lattices is sharp enough, a sub- 
stantial range of "critical" behavior of Ds — and of re- 
laxation times — might be observed near to an apparent 
transition at pa ~ Pip- This reasoning can be general- 
ized to other KA models — indeed a sharper crossover 
should occur in higher dimensions — and could be an 
explanation of the apparent dynamical transition found 
for the original cubic lattice KA ||. 

We have shown how simple kinetically constrained 
models exhibit concretely some of the qualitative 
features conjectured to be the cause of the dramatic 
slowing down in structural glasses. In particular, 
finite-dimensional models can "almost" have a a dy- 
namical phase transition; such exists on Bethe lattices 
and percolation arguments suggest it is replaced by a 
crossover — possibly very sharp — in finite dimensions. 
The dominant dynamical processes involve cooperative 
rearrangements of regions whose size diverges in the high 
density limit; nevertheless, these yield a non- vanishing 
self-diffusion coefficient at all densities. But the length 
and time scales of these processes grow extremely 
rapidly at high densities. This is reminiscent of the 
super- Arrhenius increase of the relaxation time as fragile 
liquids are cooled. 
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